Astronomy & Astrophysics manuscript no. 13860 


©ESO 2010 


January 29, 2010 





Magnetic Doppler imaging of or 2 Canum Venaticorum in 
all four Stokes parameters* 

Unveiling the hidden complexity of stellar magnetic fields 

O. Kochukhov 1 and G. A. Wade 2 

1 Department of Physics and Astronomy, Uppsala University, Box 516, Uppsala SE-75 1 20, Sweden 

2 Department of Physics, Royal Military College of Canada, Box 17000, Kingston, Ontario, K7K 4B4 Canada 

Received 14 December 2009 / Accepted 28 January 2010 

ABSTRACT 

Context. Strong organized magnetic fields have been studied in the upper main sequence chemically peculiar stars for more than half 
a century. However, only recently have observational methods and numerical techniques become sufficiently mature to allow us to 
record and interpret high-resolution four Stokes parameter spectra, leading to the first assumption-free magnetic field models of these 
stars. 

Aims. Here we present a detailed magnetic Doppler imaging analysis of the spectropolarimetric observations of the prototypical 
magnetic Ap star a 2 CVn. This is the second star for which the magnetic field topology and horizontal chemical abundance inhomo- 
geneities have been inferred directly from phase-resolved observations of line profiles in all four Stokes parameters, free from the 
traditional assumption of a low-order multipolar field geometry. 

Methods. We interpret the rotational modulation of the circular and linear polarization profiles of the strong Fe n and Cr n lines in the 
spectra of a 2 CVn recorded with the MuSiCoS spectropolarimeter. The surface abundance distributions of the two chemical elements 
and a full vector map of the stellar magnetic field are reconstructed in a self-consistent inversion using our state-of-the-art magnetic 
Doppler imaging code InversIO. 

Results. We succeeded in reproducing most of the details of the available spectropolarimetric observations of a 2 CVn with a mag- 
netic map which combines a global dipolar-like field topology with localized spots of higher field intensity. We demonstrate that these 
small-scale magnetic structures are inevitably required to fit the linear polarization spectra; however, their presence cannot be inferred 
from the Stokes / and V observations alone. We also found high-contrast surface distributions of Fe and Cr, with both elements show- 
ing abundance minima in the region of weaker and topologically simpler magnetic field. 

Conclusions. Our magnetic Doppler imaging analysis of a 2 CVn and previous results for 53 Cam support the view that the upper 
main sequence stars can harbour fairly complex surface magnetic fields which resemble oblique dipoles only at the largest spatial 
scales. Spectra in all four Stokes parameters are absolutely essential to unveil and meaningfully characterize this field complexity 
in Ap stars. We therefore suggest that understanding magnetism of stars in other parts of the H-R diagram is similarly incomplete 
without investigation of their linear polarization spectra. 
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1. Introduction 

Magnetic fields play a fundamental role in the physics of the at- 
mospheres of a significant fraction of stars on the H-R diagram. 
The magnetic fields of intermediate-mass upper main sequence 
stars (the Ap stars) have very different characteristics, and prob- 
ably a di fferent ori gin, than those of late- type stars like the 
sun (e. g. lMestefll2003h . In Ap stars, the large-scale surface mag- 
netic field is observed to be static on timescales of at least many 
decades, and appears to be "frozen" into a rigidly rotating atmo- 
sphere. The magnetic field is globally organized, permeating the 
entire stellar surface, with a relatively high field strength (typi- 
cally of a few hundreds up to a few tens of thousands of gauss). 
Magnetic field appears to be present in only a small fraction of 
intermediate-mass main sequence stars, and its presence strongly 
influences energy and mass transport (e.g., diffusion, convection 
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and weak stellar winds) within the atmosphere. The character- 
istic consequence of this interaction is the presence of strong 
chemical abundance non-uniformities in photospheric layers. 

The weight of opinion holds that the magnetic fields of up- 
per main sequence stars are primarily fossil fields - the slowly- 
decaying remnants of magnetic field accumulated or generated 
during the complex process of star formation. In this framework, 
the global and local properties of the magnetic field (the mean 
intensity, obliquity relative to the stellar rotation axis, and large- 
scale topology; but also its smaller-scale structure) may well pro- 
vide unique information about invisible interior processes which 
have occurred or are occurring within the star: differential ro- 
tation, meridional circulation currents, global- and local-scale 
dynamo action, etc. Therefore a knowledge of the detailed mag- 
netic field structure of Ap stars - on both large and small scales - 
can contribute significantly to understanding the physics of such 
stars. 

The characteristics of magnetic fields in upper main se- 
quence stars are inferred from the influence of the Zeeman ef- 
fect on their spectra. Even a relatively weak magnetic field 
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(a few tens of gauss) will clearly imprint its presence on the 
local emergent spectrum of a star, by splitting and polariz- 
ing line profiles. The large majority of magnetic field data 
in the literature measure the longitudinal Zeeman effect via 
induced circular polarization within photospheric absorption 
lines, as quantified by the mean longitudinal magnetic field. 
While such measurements represent a powerful tool for de- 
tecting organized magnetic fields, they are relatively insensi- 
tive to the topology of the magnetic field, typically constrain- 
ing the strength and orientation of the field's dipole compo- 
nent. If line circular polarization (the V Stokes parameter) is 
measured using high spectral resolving power, additional detail 
can be recovered by exploiting the rota tional Doppler effect via 
Magnetic Doppler Imaging (M DI; see Pis kunov & Kochukhovl 
2002tlKochukhov & Piskunovl2002l) . This was demonstrated by 
Kochukhov et al. (120021) . who mapped the magnetic field of the 
Ap star a 2 CVn using a timeseries of high-resolution Stokes / 
and V line profiles, interpreted using MDI. To converge to a 
unique solution, the map of a 2 CVn was forced to resemble, 
to the greatest extent possible given the constraints imposed by 
the data, a non-axisymmetric multipolar configuration. The re- 
sulting map showed a dominant dipole component, with a small 
quadrupolar contribution. Overall, the map was quite smooth, in 
agreement with expectations. 

Additional magnetic field structural detail is accessible 
through the spectral line linear polarization (Q and U Stokes pa- 
rameters) induced by th e transverse Zeeman effect. As demon- 
strated bv lLerovl Jl9 95). Zeeman linear polarization, even when 
measured photometrically using broad bandpasses, provides in- 
formation about the smaller-scale structure of th e magnetic 
field t hat is not available from circular polarization. Wad e et al.l 
(2000b) extended these observational investigations to their nat- 
ural culmination by obtaining high-precision measurements of 
Zeeman polarization in both circular and linear polarization (i.e. 
in all four Stokes /, V, Q and U parameters) for a sample of 
about 15 Ap stars, with a resolving power sufficient to distin- 
guish the variation of the polarization across individual spectral 
lines. For about a dozen of these stars, multiple Stokes IQUV se- 
quences were acquired, allowing a characterisation of the longi- 
tudinal and transverse components of the magnetic field on dif- 
ferent regions of th e st ellar surface. For the A2p star 53 Cam, 
IWade et al.l d2000bl) and[Bagnulo etaD (120011) demonstrated that 
the existing low-order multipolar models of the magnetic field 
were unable to reproduce the observed Stokes Q and U spectra, 
although they acceptably matched the Stokes V spectra. They in- 
terpreted their results as indicating that the magnetic topology of 
53 Cam was significantly more complex than that envisioned by 
the models. 

To test this proposal, Kochuk fiov et all (l2004b[) applied MDI 
to the 53 Cam data of IWade et al.l d2000bl) and Baenulo et all 
(2001) - the first attempt to reconstruct a stellar magnetic field 
using measurements of spectral lines in the four Stokes param- 
eters. Kochukhov et al. were able to reproduce the variable in- 
tensity and morphology of the polarized line profiles of 3 strong 
Fe n lines throughout t he star's ro t ational cycle. In agreement 
with t he prediction of Wad e et al.l d2000b) and iBagnulo etaD 
(2001), the derived field topology was surprisingly complex. 
In particular, while the radial field exhibited an approximately 
dipolar behaviour, the tangential field and the field modulus 
showed structure on much smaller scales. 

In the context of their analysis of 53 Cam, Kochukhov et al. 
demonstrated that the constraint on the transverse field provided 
by linear polarization measurements was essential to detecting 
the fine structure of the field. It was therefore natural to revisit 



the analysis of a 2 CVn, to examine to what extent the smooth 
map derived by iKochukhov et al.l (120021) is attributable to their 
lack of linear polarization data. This is the topic of the present 
paper. In Sect. [2] we describe the spectropolarimetric observa- 
tions from which the MDI maps are derived. In Sect. [3] we de- 
termine the physical parameters of the star, with the particular 
goal of constraining the orientation of the stellar rotation axis 
with respect to the line-of-sight. In Sect. |4]we review basic prin- 
ciples of MDI, and describe the spectral lines used, their asso- 
ciated atomic data, and the determination of global parameters 
to be used in the mapping procedure. In Sect. [5]we describe the 
results of the mapping, paying particular attention to the signif- 
icance and robustness of small-scale structures detected in the 
field modulus map. 

2. Spectropolarimetric observations 

The main set of spectropolarimetric observations of a 2 CVn 
was obtained in 1997-1999 using the (now decommissioned) 
MuSiCoS spectropolarimeter at the 2-m Bernard Lyot Telescope 
(TBL) at Pic du Midi observatory. 

The MuSiCoS instrument consists of a table-top c ross- 
dispersed echelle spectrograph (Baud rand & Bohmlll992l) . fed 
by optical fibres directly from a Cassegrain-mounted polariza- 
tion analysis module. This instrument allows the acquisition of 
a stellar spectrum in a given polarization state (Stokes V, Q or 
U) throughout the spectral range 450 to 660 nm in a single expo- 
sure. The resolving power is A/AA ^ 35 000. The optical charac- 
teristics of the spectropolarimeter an d corresponding observing 
procedures are described in detail bv lDonati et alj d!999l) . 

When MuSiCoS is used for polarization observations, 
starlight enters the polarimeter at the Cassegrain focus. The 
beam then may optionally pass through a rotatable A/ 4 re- 
tarder (in the case of Stokes V observations) or not (in the 
case of Stokes Q or U observations). The beam then inter- 
sects a Savart-type beamsplitter which separates the stellar light 
into two beams which are respectively polarized along and per- 
pendicular to the instrumental reference azimuth. The analyzed 
beams are then injected into the double 50 //m fibre, which trans- 
ports the light to the spectrograph. Spectra in both orthogonal 
polarizations are thereby recorded simultaneously on the thinned 
1024 x 1024 pixel site CCD detector. 

A single polarimetric observation, yielding an intensity spec- 
trum and one other Stokes parameter, consists of a sequence of 4 
subexposures, between which the retarder (for circular polariza- 
tion Stokes V) or the polarimetric module itself (for linear polar- 
izations Stokes Q an d U) is rotat e d by ±90° following the pro- 
cedure suggested by ISemel et al.l d!993l) . This has the effect of 
exchanging the beams within the whole instrument and switch- 
ing the positions of the two orthogonally polarized spectra on the 
CCD. This observing procedure suppresses all first-order spuri- 
ous polarization signatures down to well below the noise level. 

Spectra obtained at the TBL using the MuSiCoS spectro- 
graph and polarimeter were reduced using the ESpRIT reduction 
package (Don atiet al.ll 19971) . The ESpRIT-reduced spectra were 
then post-processed to improve the quality of the continuum nor- 
malization. A full description of the acquisition and reduction 
pr ocedure of the MuS iCoS observations of a 2 CVn is provided 
bv lWade et alJ d2000bl) . 

A complete list of the spectropolarimetric observations of 
a 2 CVn analyzed in our study is presented in Table [1] This table 
provides information on the UT date of observation, lists Stokes 
parameters obtained, their Julian dates, mean rotation phases and 
the signal-to-noise ratio. The Stokes V spectrum recorded on 25 
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February 1997 was treated separately from the linear polariza- 
tion observations on that night due to a gap in observing time 
exceeding 1 % of the rotation period. 

A close examination of the Stokes Q and U profiles ob tained 
on 20 February 1997 Qp = 0.391, see Wad e et alj|2000bl) . com- 
paring with the observations closest in phase, suggests that their 
signs are inverted. These particular linear polarization spectra 
were obtained with an interruption in the observing procedure, 
which possibly led to a sign error. We have obtained an ad- 
ditional observation of a 2 CVn on 26 July 2004, at the rota- 
tion phase <p = 0.413. The Stokes Q and U profiles from this 
phase are inverted in comparison to the Feb 20, 1997 observa- 
tion and agree with the other spectra obtained in 1997-1999. 
Consequently, we decided to disregard the linear polarization 
spectra from 20 February 1997, using only the Stokes V data 
from that night. 

Thus, in the magnetic DI analysis presented below we have 
modeled 20 separate phases of partial or complete Stokes param- 
eter observations of a 2 CVn. For the rotation phases where 2 or 
3 intensity profiles were available with each of the other Stokes 
parameters, we have constructed an average Stokes / spectrum. 
All spectra were then interpolated onto a common wavelength 
grid as required for comparison with theoretical calculations. 

All observations of a 2 CVn are phased using the ephemeris 
of|F arnswo rthl (fT93l : 

7D(£ Mmax ) = 2419869.720 + 5 d 46939 x E, (1) 

which gives the time of the maximum intensity of the spectral 
lines of Eu n and roughly corresponds to the negative extremum 
of the longitudinal magnetic field phase curve dWade et alj 
I2000ah . 

3. Physical properties of a 2 CVn 

Our modeling of the Stokes IQUV spectra of a 2 CVn and the in- 
terpretation of their rotational variability in terms of the surface 
maps of both magnetic field and chemical elements requires de- 
termining certain basic physical and geometrical properties of 
the target star: 

i Atmospheric parameters - T e g, log g and average chemical 
composition - are needed as inputs for a model atmosphere 
code, which provides us with the vertical temperature- 
density stratification used for the numerical calculation of 
the local Stokes profiles and continuum intensities; 

ii Rotation period P ml and projected rotational velocity v e sin i 
are needed to establish the mutual phasing of observations 
and to evaluate the rotational Doppler shift across the stellar 
surface; 

iii Two angles, i and 0, establish the orientation of the 
stellar axis of rotation wit h respect to the observer 
dPiskunov & Kochukhovl2002l) . 

Some of these parameters can be adjusted by optimiz- 
ing the fit to the observ ations with multiple inversions 
(iKochukhov & Piskunovll2002l) . Other parameters produce little 
impact on the Stokes spectra and Doppler images and hence re- 
quire a prior estimation. 

3.1. Atmospheric parameters 

o- 2 CVn (12 CVn A, HD 112413, HR4915, HIP 63125) 
was the first star given the "A-type peculiar" classification 



(iMau rv & Pickering] ! 18971) and is a prototype of s pectrum vari- 
able s tars. It was classified as AOpSiEuHg by ICowlev et alJ 
(1 19691) and was targeted by many notable magn e tic and spec- 
trum variability stu dies (Ba bcock & Burdl fl952t iPvperi 1969: 
ICoherj|1970t lB orra & Landstreetl 19771) . including early attempts 
to reconstruct maps of abundance s pots and magnetic field 
with the Doppler Imaging t echnique (iGoncharskii et "aD 1 1983t 
iKhokhlova & Pavlovdll984T:lGlagolevskii et al.lll985l) . 

In the previous DI study of a 2 CVn dKochukhov et alJl2002h 
we have determined the effective temperature and surface gravity 
using theoretical model atmosphere fit to the optical spectropho- 
tometry and hydrogen B aimer line profiles. We found that the 
ATLAS9 models computed with 10 times solar metal abundance 
provide a satisfactory fit to observ ations for T e ff = 1 1600 + 150 K 
and logg = 3. 9 ±0.1. Subsequentlv lLipski & Stepieiil 12008) used 
the same model atmosphere code to fit the optical spectral en- 
ergy distribution simultaneously with the UV fluxes recorded 
by the IUE satellite. They showed that the UV data require 
Teff ~ 10750 K, which is about 1000 K cooler than inferred 
from the optical spectro photometry. The f i nal T p g =11250 + 
500 K recommended by iLipski & Stepieril d2008l) is formally 
not dif ferent from the temperature adopted by IKochukhov et all 
(2002), but has a much larger uncertainty reflecting the difficulty 
of fitting the stellar spectrum in a broad wavelength region. 

It is evident that further progress in deriving a realistic 
model structure of the atmosphere of a 2 CVn can be achieved 
by computing mod els with individualised chemical compositio n 
and magnetic field dShulvak et al.l2004lKochukhov et al.l2005l) . 
Nevertheless, for the purpose of our investigation, and to retain 
the possibility of a direct comparison with the results of the 
previous DI st udy, it is sufficient to u se the model atmosphere 
parameters of IKochukhov et all (|2002), which represent reason- 
ably well the mean atmospheric properties of the star. Using for- 
ward spectrum synthesis calculations we have verified that mod- 
ification of r e ff by +500 K produces only a marginal change of 
the line intensities, comparable or smaller than the uncertainty 
due to imperfect atomic data. Thus, the temperature uncertainty 
is not a critical factor in our analysis of a 2 CVn. 

In addition to the uncertainties mentioned above, the as- 
sumption of a constant model structure is not fully correct for 
spotted Ap stars since the line and continuous opacity is differ- 
ent inside and outside spots. However, this common assumption 
is justified for Doppler imaging studies similar to ours because 
the local line profiles are sensitive to model structure effects to a 
much smaller degree than to changes of abundance or magnetic 
field. 

Taking into account the distance to a 2 CVn, d = 35.2 + 
1.1 pc which follo ws from the revised Hipparcos pa rallax 
(Ivan Leeuw en 2007), and adopting V = 2.89 + 0.04 from Pvper 
( 1969), we find the absolute visual magnitude M v = 0.16 + 6.08. 
For the temperature T e g= 11600 + 5 00 K and the bolomet- 
ric correction BC = -0.42 measured bv lLipski & Stepienl d2008l) 
we obtain L- 101 + 12 L H and R = 2 A9 ± 0.26 R Q assuming 
M® o] =4.75 for the Sun dBessellll2000h and using an uncertainty 
of 0. 1 for the bolometric correction of a 2 CVn. Then, employ- 
ing this stellar luminosity and temperature we can determine the 
evolutionary state of a 2 CVn by comparing its observed posi- 
tion in the H-R diagra m with the prediction s of the theoretical 
evolut ionary tracks by [Schaller et al.1 d 19921) and ISchaerer et alj 
d!993h . Assuming an overall metallicity Z = 0.018, we derive 
a stellar mass M = 2.97 + 0.07 M and an age of 1.65+°!; x 10 8 
years. These parameters yield the surface gravity log g = 4.12 ± 
0.09, which can be brought in agreement with the spectroscopic 
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Table 1. Four Stokes parameter spectropolarimetric observations of a 2 CVn used for Magnetic Doppler Imaging. 



UT Date 




Stokes 




JD (2450000+) 






S/N 




Parameters 


Q 


U 


V 


f 


6<p 


range 


18 Feb 1997 


Q 


U 


V 


498.511 


498.531 


498.489 


0.038 


0.004 


420-550 


19 Feb 1997 


Q 


u 


V 


499.532 


499.555 


499.511 


0.225 


0.004 


650-690 


20 Feb 1997 a 






V 






500.422 


0.387 




550 


21 Feb 1997 


Q 


u 


V 


501.490 


501.504 


501.470 


0.582 


0.003 


270-320 


22 Feb 1997 


Q 




V 


502.549 




502.519 


0.773 


0.003 


320-490 


23 Feb 1997 


Q 


u 


V 


503.497 


503.516 


503.474 


0.949 


0.004 


240-280 


25 Feb 1997 b 






V 






505.581 


0.331 




470 


25 Feb 1997 


Q 


u 




505.695 


505.715 




0.353 


0.002 


410-410 


07 Feb 1998 


Q 


u 


V 


852.562 


852.584 


852.539 


0.771 


0.004 


740-890 


09 Feb 1998 


Q 


u 


V 


854.561 


854.584 


854.539 


0.137 


0.004 


820-1020 


11 Feb 1998 


Q 


u 


V 


856.545 


856.567 


856.524 


0.499 


0.004 


910-1020 


12 Feb 1998 


Q 


u 


V 


857.675 


857.694 


857.657 


0.706 


0.003 


820-880 


lj rcu I770 


Q 


u 


V 


858.555 


858.574 


858.537 


U.oO / 




7in sin 

/ 1 U— o 1 u 


15 Feb 1998 


Q 


u 


V 


860.605 


860.626 


860.586 


0.242 


0.004 


500-720 


16 Feb 1998 


Q 


u 


V 


861.606 


861.627 


861.584 


0.424 


0.004 


560-750 


17 Feb 1998 


Q 


u 


V 


862.533 


862.554 


862.515 


0.594 


0.004 


560-560 


19 Feb 1998 


Q 


u 


V 


864.603 


864.638 


864.576 


0.973 


0.006 


430-590 


15 Jan 1999 


Q 




V 


1194.621 




1194.600 


0.310 


0.002 


400-450 


19 Jan 1999 


Q 


u 


V 


1198.652 


1198.672 


1198.633 


0.049 


0.004 


770-840 


26 Jul 2004 


Q 


u 


V 


3213.362 


3213.381 


3213.400 


0.413 


0.004 


440-480 



Columns give the UT date for the beginning of the night when observation was obtained, the Stokes parameters observed (in addition to Stokes 
/), the Julian Dates of mid-exposure for each observed Stokes parameter, the mean phase 7p for all observed Stokes parameters, the maximum 
difference Sip between 7p and the phases of individual Stokes parameter observations, the range of the signal-to-noise ratio per 2.6 km/s spectral 
pixel for two or three exposures. 

a The Stokes Q and U spectra are disregarded. See text for detailed explanation. 

b The Stokes V observation is considered individually due to a larger-than-average phase difference with respect to the linear polarization spectra 
obtained on that night. 



Table 2. Fundamental parameters of a 2 CVn. 



Parameter 



Value 



Ref. 



7eff 


11600 + 500 K 


(1), (2) 


l°gg 


3.9 + 0.1 


(1) 


Distance 


35.2+1.1 pc 


(3) 


My 


0.16 + 0.08 


(4) 


BC 


-0.42 + 0.1 


(2) 


L 


101 ± 12 L 


(4) 


R 


2.49 ± 0.26 R e 


(4) 


M 


2.97 ± 0.07 M Q 


(4) 


Age 


1.65 + ° * x 10 s yr 


(4) 


v e sin i 


18.4 ± 0.5 km s- 1 


(4) 


*rot 


5?46939 


(5) 


i 


120 ± 5° 


(4) 


& 


115 + 5° 


(4) 



Re ferences. (l)[K ochukho v et alj i2002h : ( 2)lLipski & Stepienl d2008l) : 
(3) lvan Leeuwenl j2007h ; (4) this work; (5) lFarnsworthUl932h . 



estimate if the latter is corr ected for the stro ng He deficiency of 
the atmosphere of a 2 CVn dKochukhov et al.ll2002h . 

The fundamental parameters of a 2 CVn determined in our 
study, as well as those adopted from previous investigations, are 
summarized in Table [2] 

3.2. Rotational parameters 

The widely- used rotation perio d of a 2 CVn, 5^46939, was de- 
termined by Farnsworth] (Il932l) based on a series of intensity 
measurements of the blue Eu n lines taken at several observato- 
ries during the years 1912-1932. Using this rotational period one 
can successfully phase together longitudinal magnetic field mea- 



surements obtained over the time spa n of 25 year s ( Wade et ail 
2000a]). Thus, the rotational period of Farnsworth] (1 1932b is ade- 
quate for our analysis. We do not attempt to improve it since no 
recent extensive photometric or spectroscopic time-series obser- 
vations are available for a 2 CVn. 

The projected rotational velocity is one of the key input pa- 
rameters for DI inversio ns and needs to be determined with a rel- 
atively high precision dKochukhov & Piskunovl l2002). It is best 
found by optimizing the fit to line profiles observed at high spec- 
tral resolution. As described in Sect. 14.31 the intensity profiles 
of the strong Fe n and Cr n lines in the spectrum of a 2 CVn are 
reproduced with v e sin i = 18.4 + 0.5 kms -1 . T his value is some- 
what l arger than 17.4 + 0.5 kms -1 found bv iKochukhov et all 
(2002) using Cr n lines of different strength. Although formally 
not significant, this discrepancy could reflect the difference in in- 
strumental profile or, possibly, be an artifact of vertical chemical 
stratification neglected in our spectrum synthesis. 

The orientation of the stellar rotational axis is characterized 
by the two angles i and 0. The inclination, 0° < i < 180°, is 
the angle between the rotational axis and the observer's line- 
of-sight. The values i > 90° correspond to the situation when 
we see a clockwise rotation of the star from the visible ro- 
tation pole. The azimuth angle, 0° < < 360° determines 
the sky-projected position angle of the rotational axis. This an- 
gle is counted coun terclockwise from the North Celestial Pole 
(lLandolfi et al.|[l993h . The azimuth angle is relevant and the dis- 
tinction between i and 180° - i can be made only when interpret- 
ing linear polarization observations. Note that since the Stokes 
Q and U parameters depend on trigonometric functions of 20, 
the values of and © + 180° cannot be discriminated. 

Num erical tests of the four Stokes parameter inversion pro- 
cedure dKochukhov & Pi skunov 2002) indicate that the tilt and, 
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Fig. 1. Photopolarimetric measurements of the longi t udinal mag- 
netic field in a 2 CVn dBorra & Landstreed 119771: lLandstreetl 
1982, upper panel) and the total net linear polarizatio n derived 
from th e LSD Stokes Q and U measurements of Wad e et al.l 
( 2000a) (lower panel) compared with the predictions of the dipo- 
lar model (solid curves) with parameters B p = 4.6 kG, i = 119° 
and B = 78°. The model with the same B p and i' = 180° - i, 
B' = 180° - B yields indistinguishable curves. 



especially, the azimuth angle of the rotational axis can be re- 
liably determined by minimizing the x 2 of the fit to observa- 
tions in inversions with different i and 0. However, our magnetic 
DI technique is fairly costly in computing time, which makes 
exploring all possible rotational axis orientations impractical. 
Instead we start by finding a plausible range for each angle us- 
ing a simple dipolar field model and then refine the angles with 
a series of MDI reconstructions. 

The inclination angle is inferred from the simultaneous 
fit of the dipolar field model to the Ba lmer line photopo- 
larimetric longitudinal field measurements (Bor ra & Landstreel 
119771; lLandstreetlll982h and the total linear polarization, P L = 

^1(0.1 1) 2 + (U/I) 2 , obtained from th e net linear polarization 
measurements of IWade et al.l d2000al) . We employ longitudinal 
field measurements obtained from H lines to avoid the strong 
modulation due to to abundance nonuniformities that is evident 
in longitu dinal field measur ements derived from, e.g., lines of 
Fe or Cr dWade et al.ll2000al) . The dipole model parameters in- 
clude the polar field strength B p , magnetic obliqui ty B and angles 
i and . The theoretical P^ is computed following Landolfi et al. 
( 1993). Using the total linear polarization instead of individual 
net Q and U measurements allows us to avoid solving for at 
this stage of the analysis. 

The comparison of the magnetic observables with the dipo- 
lar model predictions is presented in Fig. [TJ The scatter around 
the best-fit curve is considerable but the model evidently repro- 
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Fig. 2. Comparison of the net linear polarization (symbols) ob- 
tained from the LSD Stokes Q and U profiles of a 2 CVn (upper 
and middle panels, respectively) with the dipolar model predic- 
tions (lines). The solid line corresponds to the magnetic field ge- 
ometry given by B v = 4.6 kG, i = 1 19°, B = 78°, = 1 10°. The 
dashed line shows predictions of the best-fitting model for the 
same B v , 180° -z, 180° -B and = 1 19°. The lower panel shows 
the total chi-square of the fit to net linear polarization curves as 
a function of the azimuth angle for the dipolar model with 
i = 119°, B = 78° (solid line) and for the model with comple- 
mentary angles 180° - i and 180° -B (dashed line). 



duces the main features of the (B z ) and P^ phase dependence if 
we adopt parameters B p = 4.6 ± 0.5 kG, B = 78 ±7°, i = 1 19 ± 12° 
(or 180° -y8=102° and 180° - / = 61°). The inclination angle 
obtained in this way is consistent with i = 53 ± 8° that can be de- 
rived from P rot , R and v e sin i of a 2 CVn using the usual oblique 
rotator relation. 

In the second step we find the values of which provide 
the best match to the observed rotational variation of the Q and 
U net linear polarization for the previously determined dipolar 
model parameters. Fig.|2]demonstrates that a global^- 2 minimum 
is found for * 1 10° and that B = 78°, i = 1 19° provide a better 
match to the observations t han 180° - B and 18 0° - i. The same 
conclusion was reached bv IWade et al.l d2000al) . Thus, prelimi- 
nary modeling suggests i w 120° and » 110°. The final values 
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of these parameters are established in Sect. I4.3l with the help of 
multiple MDI inversions which take into account chemical in- 
homogeneities and incorporate deviations of the magnetic field 
topology from a simple dipole. 

4. Magnetic Doppler imaging 

4.1. Magnetic inversion technique 

A detailed discussion of Magnetic Doppler Imaging using high- 
res olution circular and linear pol arization spectra is provided 
bv lPiskunov & Kochukhov! (I2002 I). In that publication the MDI 
code InversIO was introduced, and key numerical techniques 
employed in magnetic inversions were discussed. Here we 
briefly recount the b asic principles of our MD I meth odology, re- 
ferring the reader to iPiskunov & Kochukhov! d2002l) for further 
details and for an in-depth description of technical issues. 

The aim of the magnetic inversion is to derive the surface 
distributions of magnetic field and chemical abundances by min- 
imizing the total discrepancy function 



¥ = D + % 



(2) 



where T) characterizes the discrepancy between the observed 
and computed phase-resolved spectra and % is the regularization 
functional. For DI with polarization data 



•■) 



rrrobsl 2 / 2 



(3) 



where Tjgpi are the computed and observed Stokes parameter 
profiles for rotation phase <p and wavelength point A. The in- 
dex k represents summation over the available Stokes param- 
eters (all four in the present study of or 2 CVn), while weights 
ojic are introduced to ensure approximately the same contribu- 
tion of each Stokes parameter to D. In practice this is achieved 
by choosing <x>k proportional to the phase-averaged amplitudes 
of the observed Stokes parameters. For our a 2 CVn data set 
ll>i : (Dq : u>u : «y = 1 : 45 : 45 : 6. 

Theoretical Stokes IQUV spectra for a given rotation phase 
depend on the surface topology of the magnetic field B and the 
geometry of the abundance distributions s , s 2 , etc. Normalized 
Stokes profiles are obtained by summing the Doppler-shifted lo- 
cal Stokes spectra X\ over a discrete surface grid and dividing 
by the disk-integrated intensity in the unpolarized continuum 



(4) 



Here, the index i runs over a set of N v surface zones visible to the 
observer at the rotational phase ip and S- is the projected area of 
surface element i. For the present study of a 2 CVn we use a grid 
with 695 surface zones, which is fully sufficient given the mod- 
erate v e sin i of the star and the limited spectral resolution of the 
observational data. The local Stokes parameter profiles Xf* are 
computed for a given model atmosphere and local values of the 
abundances and magnetic field by solving numerically the polar- 
iz ed radiative transfer equation in its general form, as described 
bv lPiskunov & Kochukhov! J2002). Thus, theoretical spectra are 
computed by InversIO fully self-consistently with the current 
magnetic and abundance maps and treating the line formation in 
the magnetized stellar atmosphere without any of the simplifying 
assumptions that are commonly used, e.g. a Milne-Eddington at- 
mosphere or fixed local Gaussian line profiles. 



We apply regularization to ensure stability of the complex 
optimization process of magnetic DI and to obtain a unique so- 
lution independent of the initial guess and of the surface dis- 
cretization. The general form of the regularization functional % 
implemented in InversIO is 



A i Z Z (B < - b j )2 + Ai Z (Bi - ^ ult)2 



(5) 



+ A, 



ZZ^-*)) 2 +ZZ^ 2+ 



The first and third terms represent the Tikhonov regularization 
functions for magnetic field and abundance distributions, re- 
spectively. Ai and A3 are corresponding regularization param- 
eters, selected in such a way that DfR^ 2-10 at the conver- 
gence. Application of the Tikhonov regularization leads to lo- 
cal smoothing of the solution by minimizing the difference be- 
tween the value of parameters in the grid point i and all neigh- 
bouring points /'. Appropriate choice of the Tikhonov regulariza- 
tion ensures that the magnetic and abundance distributions re- 
constructed by our code represent the simplest possible solution 
which is still compatible with observational data. At the same 
time, the Tikhonov regularization does not impose an a priori 
mathematical model for the global magnetic field geometry or 
abundance distribution. 

A comprehensive ser ies of numerical exper i ments with 
InversIO was presented by Kochu khov & Piskunovl (l2002). This 
paper explored the performance of the magnetic DI method 
for different combinations of stellar parameters, magnetic field 
strengths and geometries by reconstructing surface maps from 
simulated Stokes parameter observations. These numerical tests 
convincingly demonstrate that a successful inversion, reveal- 
ing magnetic field structures at both the large and small scales, 
is possible when the code is applied to the full Stokes vec- 
tor d ata set in conju nction with Tikhonov regularization. In 
fact, iPiskunovl d2005l) showed that, from the basic mathemati- 
cal standpoint, the underlying inverse problem of magnetic map- 
ping with all four Stokes parameters is well-posed and has a 
unique solution. In this case regularization is only needed to 
suppress instabilities produced by the noise in real data and its 
sparse phase and wavelength coverage. However, due to the lack 
of high-quality Stokes IQUV stellar spectra, the only previous 
magnetic DI inversion based on polarization observations in all 
four Stokes parameters was our study of the Ap star 53 Cam 
dKochukhov et alj|2004al) . 

While Stokes IQUV datasets suitable for Magnetic DI are 
relatively rare, comparable Stokes IV datasets are more com- 
mon. Magnetic DI with only circular polarization observations 
has therefore been more widely applied. However, such datasets 
result in maps which are less informative and robust com- 
pared to those computed from full Stokes vector data. The IV 
mapping problem is intrinsically ill-conditioned and the out- 
come of such inversions depends sensitively on the choice 
of regularization function and the initial guess. Our numeri- 
cal experiments with InversIO dPiskunov & Kochukhov 2002; 
iKoch ukhov & Pisk unovl 120021) showed that a successful re- 
construction of the global field topology of Ap stars using 
only Stokes IV data requires imposing a prior assumption that 
strongly restricts the range of possible derived magnetic field 
topologies. In our DI code this assumption is implemented in 
the form of so-called multipolar regularization (the second term 
in Eq. ©), which directs the solution toward a general, non- 
axisymmetric second-order multipolar expansion of the mag- 
netic field structure. The multipolar model field B mult , calculated 
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Table 3. Atomic parameters of spectral lines employed in mag- 
netic inversions. 
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-0.970 


4.5 


4.5 


1.34 


1.34 


Fen 


4824.198 


10.288 


-1.215 
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4.5 


1.45 


1.22 


Fen 


4824.838 


8.145 


-1.893 
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1.25 


1.24 


Fen 
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2.891 


-1.320 
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2.00 


2.40 


Fen 


5018.440" 


2.891 


-1.271" 
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2.00 


1.87 


Fen 


5019.462 


5.569 


-2.697 


3.5 


4.5 


1.14 


1.10 



The columns specify the ion, the central wavelength A, excitation 

potential of the lower atomic level E\„, quantum numbers J and Lande 

factors g for both the lower and upper levels. 

a Dominant blend components. 

b Oscillator strength is adjusted in the MDI inversion. 



by our code at each iteration, is equivalent to the dipole plus 
quadrupole parameterization of Bagnul oet al.l(ll996l). Stokes IV 
time series of the Ap stars a 2 C Vn (jKochukh ov et al. 20021). 
HP 2 47 12 dLuftinger et alJl2010h and HP 72106 dFolsom etafl 
2008) were interpreted using InversIO in the multipolar regu- 
larization mode. All these studies inferred basically dipolar field 
geometries for the target stars. 

The availability of MuSiCoS four Stokes parameter time- 
resolved spectra for a 2 CVn provides us with the capability to 
perform magnetic inversion constrained only by local, Tikhonov 
regularization (Ai + 0, A2 = 0, see Sect. 15. II ). At the same time, 
it is also very instructive to compare these ultimate MPI results 
with the maps obtained disregarding linear polarization obser- 
vations and globally constraining the magnetic inversion with 
the multipolar regularization. In Sect. 15. 21 we present results for 
magnetic PI in the circular polarization mode, with Ai =0, 
A 2 * 0. 



4.2. The choice of spectral lines 

The amplitude of linear polarization signatures in metal lines of 
even the most s trongly-magnetic Ap stars is typically at the level 
of only 10~ 3 (Wade et al. 2000b). Consequently, we are limited 
to studying the Stokes Q and U profiles of a few deepest, mag- 
netically sensitive individual lines in the stellar spectra. Using 
such features is usually not optimal for abundance mapping be- 
cause the lines are saturated over a significant part of the stel- 
lar surface and thus their variation is subdued in comparison to 
weaker lines of the same ions. For this reason, and also due to 
the relatively low resolving power of the MuSiCoS observations, 
we do not expect to achie ve the same quality of ab undance inver- 
sions as demonstrated by iKochukhov et alj d2002l) in the previ- 
ous PI analysis of a 2 CVn. However, this reduced sensitivity of 
the strong lines to abundance inhomogeneities make them more 
suitable for our primary goal of the mapping magnetic field us- 
ing four Stokes parameter profiles. 

We found three lines in the spectrum of a 2 CVn appropri- 
ate for magnetic PI. The two Fen features at A 4923.93 and 
5018.44 A were previously used bv IKochukhov et al.1 ([2004a) 
and lKhalack & Wadg d2006l) in the studies of other Ap stars ob- 
served with the MuSiCoS spectropolarimeter and proved to be 
effective for the reconstruction of the magnetic field topology. 
Their mean Lande factors, g — 1 .7-1 .9, are among the largest for 
metal lines of such intensity. Both lines show a clear and com- 
plex signal in the Stokes Q and U spectra, detected at high sig- 
nificance levels for the majority of rotation phases (see Fig. [5]). 
We also use the Cr 11 4824. 1 3 A line for magnetic mapping. This 
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Fig. 3. The relative standard deviation of the fit to Cr 11 and Fe 11 
Stokes / profiles is plotted as a function of v e sin i (symbols). 
The solid curve shows the quadratic fit in the 17.5-19.5 kms" 1 
region. The vertical dashed line corresponds to the optimal pro- 
jected rotational velocity v e sin i = 18.4 + 0.5 km s _1 . 



line shows less prominent Stokes Q and U profiles because it 
is somewhat weaker than the Fe n lines, and exhibits a smaller 
Zeeman effect (g = 1.3). Nevertheless, simultaneous recovery 
of the surface magnetic field distribution from the lines of more 
than one chemical element alway s improves the robustness of 
MPI (IKochukhov & Piskunovll2002l) . 

The complete line list adopted in the magnetic inversions of 
a 2 CVn is presented in Table [3] Parameters of the two strong 
Fen lines, one Cm line and thre e weaker Fe 11 ble nds are ex- 
tracted from the VALP database (Kup ka et alll 19991) . The oscil- 
lator strength of the lines from Fe 11 multiplet 42, to which the 
4923.93 and 5018.44 A lines belong, are notorious for the large 
scatter of log gf recommended in different literature sources. 
Even the relative log gf values of these lines are not known with 
the precision ne cessary for detailed spec trum synthesis mod- 
eling, prompting IKochukhov et all d2004al) to perform separate 
MPI inversions for each of these lines in their study of 53 Cam. 
To overcome this problem here we modified InversIO, adding 
the oscillator strength of the Fen 5018.44 A line to the list of 
free parameters determined by the code. Table [3] lists the final 
log gf for this lines, which is 0.051 dex lower than the value 
recommended by VALP. 

Additional line data required for calculation of the Zeeman 
splitting includes Lande factors of the upper and lower atomic 
levels and corresponding J quantum numbers. These quantities 
are also specified in Table[3] wi th Lande factor s originating from 
the theoretical computations bv lKuruc^ d 1993b and extracted us- 
ing the VALP interface. 

4.3. Optimization of v e sin i, i and © 

Magnetic inversion is sensitive to the adopted projected ro- 
tational velocity and the three-dimensional orientation of the 
stellar rotational axis. Numeri cal experiments presented by 
IKoch ukhov & Piskunov (2002) showed that solutions with in- 
correct values of these parameters tend to have a higher x 2 of 
the fit to the observed Stokes profiles. For v e sin i this degrada- 
tion of the fit quality is most clearly seen in Stokes /. A measure 
of the final discrepancy between the observed and computed in- 
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Fig. 4. The relative total discrepancy function is plotted against 
the azimuth angle adopted for the inversion. Different symbols 
show results obtained for three values of the inclination angle i. 
From these curves we infer the optimal values © = 115 + 5° and 
i = 120 ±5°. 

tensity profiles of the Cm and Fen lines is illustrated in Fig. [3] 
Examining the inversion results for 7 trial values of v e sin i in the 
range from 17 to 20 km s _I we find a minimum for v e sin i = 18- 
19 kms -1 . A parabolic fit to these results yields a minimum at 
v e sin i = 18.4 + 0.5 km s _1 . We adopt this rotational velocity for 
all subsequent inversions presented in our paper. 

We followed a similar strategy to optimize the choice of the 
inclination and azimuth angles. The modeling of the mean longi- 
tudinal field and net linear polarization suggests that these angles 
are close to 120° and 1 10°, respectively. To determine the values 
of i and @ we carried out inversions for the ranges i = 110— 
130° and = 90-140° with a 10° step for both angles. The 
normalized total discrepancy function obtained by InversIO at 
the convergence is illustrated in Fig. [4] for the entire grid of 18 
inversions. It is evident that i = 120° gives a smaller discrepancy 
than either i = 1 10° or 130°, while the optimal © is constrained 
to the 1 10-120° interval independent of inclination. Finally, we 
establish i = 120+5° and 0=1 15+5° and adopt this orientation 
of the stellar rotation axis for all magnetic inversions presented 
below. 

5. Results 

5.1. Full Stokes vector inversion 

We reconstruct the magnetic field topology and surface dis- 
tributions of Cr and Fe in a simultaneous inversion, using 20 
phases of observations in all four Stokes parameters. The best- 
fitting theoretical profiles obtained by InversIO at the conver- 
gence are shown in Fig. with the thick solid line. Evidently, 
the observed line profile variability in the intensity and polar- 
ization is satisfactorily reproduced. The only minor systematic 
discrepancy is seen in the profile of the Fen 5018.44 A line, 
which is overestimated in the theoretical spectra correspond- 
ing to the rotation phases around negative magnetic extremum 
(<p = 0.949-0.049). The linear polarization profiles are well-fitted 
by InversIO. Occasional discrepancies between the observations 
and the model spectra (e.g., Stokes Q profile of Fen 5018.44 A 
for (f = 0.225) are always limited to only one out the three lines 
studied and thus is likely due to observational uncertainties. 



Especially impressive is the ability of the model to reproduce 
the systematic details of the linear polarization profiles, such as 
the peculiar "w"-shape of Stokes U around phase 0.4. The agree- 
ment between the observations and synthetic Stokes V spectra is 
excellent. 

The surface magnetic field distribution of a 2 CVn recovered 
by InversIO from the phase-resolved Stokes profiles of the Cm 
and Fen lines is presented in Fig. [6] This plot shows spherical 
projections of the magnetic map, displaying separately the field 
strength and the field orientation distributions. In this and similar 
plots the star is shown at 5 equidistant rotational phases, at the 
aspect angle corresponding to i = 120° and = 0°. The overall 
structure of the magnetic field in a 2 CVn is dipolar-like in the 
sense that approximately half of the stellar surface is covered 
with the outward-directed radial field while another half exhibits 
inward-directed field. This global magnetic field topology agrees 
with the orientation of the dipolar axis suggested in previous 
magnetic analyses of a 2 CVn where m ultipolar models were fit- 
ted to integral m agnetic observables dBorra & Landstreell 19771: 
i Gerfhet all 1999b or to Stokes IV line profiles (iKochukhov eta?! 
120021) . However, the field strength distribution illustrated in the 
upper panel of Fig. [6] reveals that the derived magnetic field 
topology of a 2 CVn resembles that of a dipole on the largest spa- 
tial scale. The outcome of our four Stokes parameter inversion 
suggests that the field is, in fact, far more complex. In particular, 
there is a definite asymmetry in the field strength and structure 
between the negative magnetic pole {tp as 0.0) and the positive 
one (<p as 0.5). The field is clearly stronger in the latter case and 
its structure is dominated by the high-contrast magnetic spots 
where the field strength reaches 4.5 kG locally. This is approxi- 
mately 1 kG higher than the local field at the opposite side of the 
star. 

The rectangular projection of the magnetic field distribution 
reconstructed for a 2 CVn is presented in Fig.|7]for all three com- 
ponents of the field vector. This figure again gives an impression 
of the global dipolar topology locally distorted by small-scale 
features. The radial field component is most similar to the sur- 
face magnetic structure expected for an oblique dipole, while the 
meridional and azimuthal components deviate significantly from 
a simple, low-order multipolar shape. 

Is it possible that the complexity seen in the MDI maps 
of a 2 CVn is an artifact produced by the noise in observations 
combined with an insufficiently strong regularization adopted 
in the magnetic mapping? We have thoroughly examined this 
possibility by conducting inversions with different values of the 
Tikhonov regularization parameter for magnetic field, Ai . These 
tests confirm that our choice of regularization is correct because 
any substantial increase of the smoothing of the magnetic maps 
by the Tikhonov regularization noticeably worsens the agree- 
ment between theoretical spectra and the linear polarization ob- 
servations. For example, with the Tikhonov regularization in- 
creased by a factor of 10 the code converges to a markedly sim- 
pler field topology, with the field strength not exceeding 3 kG 
(Fig. [8). The resulting fit quality is unchanged for Stokes / and 
is only slightly worse for Stokes V relative to the profiles cor- 
responding to the complex field model. However, the smoothed 
field structure fails to reproduce the Stokes Q and U spectra at 
phases <p = 0.353-0.594, when the most complex part of the sur- 
face magnetic field distribution crosses the line of sight. For all 
these phases the observed linear polarization profiles exhibit a 
double-wave structure with two minima and two maxima across 
each spectral line (e.g. the "w"-shaped profiles of Stokes U). On 
the other hand, these features are well-reproduced by the best- 
fitting complex magnetic field distribution. In contrast, theoret- 
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ical profiles corresponding to the smoothed magnetic map lack 
the double-wave structure. The resulting degradation of the fit 
quality for the Q and U spectra is readily seen in the increase of 
X 2 by 30-60%. 

This analysis shows that the small-scale features in our mag- 
netic map of a 2 C Vn are directly connected to the particular mor- 
phology of the linear polarization signatures, observed consis- 
tently in all spectral lines studied. Thus, we can assert that our 
final magnetic field topology of a 2 CVn (Fig. [6} is truly the sim- 
plest possible field structure consistent with the available Stokes 
parameter observations. Therefore, the inferred complexity and 
substantial local deviation of the stellar magnetic geometry from 
the low-order multipolar field are real. It is remarkable that we 
are able to detect and characterize these complex magnetic struc- 
tures using Stokes IQUV data, but not using Stokes IV data. 

The surface abundance maps of Cr and Fe reconstructed by 
InversIO simultaneously with the magnetic field topology are 
presented in Fig. [9] We infer a high-contrast distribution for 
both elements, with the abundance varying by 3—4 dex and pro- 
nounced abundance minima at phase <p — 0. These abundance im- 
ages are qualit atively similar to t he iro n-peak element maps re- 
constructed by Koc hukhov et al.l (]2002) but have a lower surface 
resolution for the reasons discussed in Sect. 14.21 In addition, the 
map obtained from the two strong Fe n lines exhibits noticeably 
smaller abundance values in the region around the negative mag- 
netic pole (lo g jVFe/Wtot < -6 comp ared to logA^e/Mot ^ -4.5 
obtained by Kochuk hovet al.ll2002l) . This can possibly be as- 
cribed to the effects of vertical chemical stratification, which 
we do not take into account in our analysis. A vertical abun- 
dance distribution with a transition from a high abundance in 
the deeper atmosphere to a lower element concentration in 
the higher layers is often o bserved in cooler Ap stars (e.g., 
Kochukho 7etaT1l2006L l"2009>. If such stratification is present in 
some parts of the surface of a 2 CVn, it would lead to substantial 
weakening of the cores of intrinsically strong Fe n lines relative 
to intrinsically weak lines, yielding a smaller abundance if such 
lines are analyzed neglecting chemical stratification. However, 
these effects generally do not significantly influence polariza- 
tion profiles and thus a detailed treatment of vertical abundance 
inhomogeneities is beyond the scope of our study. 

5.2. Stokes IV inversion with multipolar regularization 

Magnetic inversion in all four Stokes parameters, similar to our 
analysis of a 2 CVn, remains very uncommon in the research 
field of stellar magnetism due to the considerable difficulties as- 
sociated with the acquisition of the necessary observational data. 
Instead, magnetic DI and Zeeman Doppler Imaging (ZDI) stud- 
ies employing time-resolved Stokes IV spectra are much more 
widely applied to study field topologies of magnetic Ap stars and 
active late-type stars. In this section we assess how much of the 
complexity of the magnetic field of a 2 CVn, evident from the full 
Stokes vector reconstruction, can be retrieved in the analysis of 
an incomplete, circular polarization-only data set. As described 
in Sect. 14.11 in this experiment we use InversIO in Stokes IV 
inversion mode with the multipolar regularization. We carry out 
magnetic field reconstruction for the same three Cr and Fe spec- 
tral features, although in principle the number of usable spectral 
lines is much larger if we are not concerned with the analysis of 
the Stokes Q and U spectra. 

The fit to the / and V spectra of a 2 CVn and prediction for 
the linear polarization profiles obtained by our code at the con- 
vergence of the IV mapping problem is illustrated with the thin 
line in Fig. [TO] To enable a direct comparison with the previ- 



ous results, this figure also shows the theoretical profiles for our 
optimal four Stokes parameter inversion (same as in Fig. 0. At 
the resolution of the MuSiCoS spectra we see no difference in 
the intensity and circular polarization between the two sets of 
synthetic profiles. Yet the picture is very different for Stokes Q 
and U. The model topology inferred from the circular polariza- 
tion alone evidently yields systematically higher amplitude of 
the linear polarization profiles and cannot match the complex 
shape of the observed Q and U signal around rotational phase 
(f « 0.5. This degradation of the fit quality is readily seen with a 
factor of 2-3 increase of the^- 2 for both Q and U spectra. Thus, a 
successful description of the phase variation of the intensity and 
circular polarization spectra gives no guarantee that the result- 
ing model of the stellar surface magnetic field is also adequate 
for the Stokes Q and U data. 

The magnetic map derived in the Stokes IV inversion is 
shown in Fig. QT| This field structure is very close to a dipo- 
lar topology with a mild non-axisymmetric quadrupolar contri- 
bution (B p = 3.7 kG, fi q = 1.0 kG). The global field compo- 
nent and distribution of the radial field is similar to those seen 
in Fig. [6] but there is no evidence of the small-scale magnetic 
structure which is needed to fit the linear polarization profiles. 
We conclude that the Stokes IV analysis of the stellar magnetic 
topologies is fundamentally limited with respect to the scale of 
magnetic structures which can be resolved with this method. In 
the particular case of a 2 CVn it can be used to assess the over- 
all, global dipolar-like field topology but, compared to the four 
Stokes parameter inversion, it cannot provide detailed maps of 
the distributions of the magnetic field strength and orientation 
across the stellar surface. 



5.3. Multipolar expansion of magnetic maps 

Spherical harmonic expansion provides a convenient method for 
quantitative assessment and detailed characterization of the mag- 
netic field maps obtained in the MDI inversions. It also allows 
us to objectively compare the magnetic field topology of a 2 CVn 
with the field geometry of other Ap stars, in particular 53 Cam 
which was previously stu died in all four Stokes parameters by 
iKochukhov et al.l (|2004a) using the same inversion technique. 

Here we apply the spherical harmonic e xpansion method in- 
troduced by Piskun ov & Kochukhov! d2002|) and employed in the 
analysis of 53 Cam by IKochukhov et al.l d2004al) . The three vec- 
tor components of the magnetic field distribution are approxi- 
mated with a superposition of real spherical harmonic series, in- 
cluding both poloidal and toroidal expansion terms. In the anal- 
ysis of a 2 CVn we truncate the multipolar expansion at I — 10. 
This gives us 240 poloidal and toroidal multipolar coefficients, 
which are determined by solving a linear least-squares problem. 
Different latitudes on the stellar surface are treated with differ- 
ent weights to account for the variation of the field reconstruc- 
tion quality and to exclude the invisible part of the stellar surface 
(hidden to us due to the stellar geometry). 

The results of the multipolar expansion are presented in 
Fig.[T2]for the magnetic field topology of a 2 CVn obtained from 
all four Stokes parameters (Fig. [6), from the Stokes IV spec- 
tra modeled using multipolar regularization (Fig. [TTJ and for the 
surface magnetic field of 53 Cam. For the latter star the MDI map 
was produced with the same modified InversIO code as applied 
in the present study of a 2 CVn, using the Stokes IQUV pr ofiles 
of the three Fe n lines studied by IKochukhov et akl (|20O4a). The 
updated magnetic map of 53 Cam does not differ appreciably 
from the average field topology inferred by those authors. 



10 



O. Kochukhov and G. A. Wade: Magnetic Doppler imaging of a 1 Canum Venaticorum in all four Stokes parameters 



Multipolar expansion coefficients, J (C"') 2 + (C ( "') 2 , are 

plotted in Fig. [12] as a function of t and m, separately for the 
poloidal and toroidal expansion terms. The greyscale for each 
magnetic map is renormalized so that the largest coefficient is 
black and the smallest is white. Fig. [12] suggests that the field 
geometry of a 2 CVn is dominated by £ — 1 modes. At the same 
time, the contribution of the higher-/' modes and the amplitude 
of the toroidal field components is also non-negligible for the 
MDI map obtained from all four Stokes parameters (left column 
in Fig. [T2l . For this surface magnetic field distribution we see 
contributions of modes with £ up to 5-6. On the other hand, the 
IV mapping results (middle column in Fig.fTZl) suggest a consid- 
erably simpler field, showing very clearly the reduction of the in- 
formation content when linear polarization spectra are excluded 
from the Doppler Imaging analysis. In this case the dominant 
dipolar mode is distorted only by a marginal contribution of the 
non-axisymmetric £ = 2 and axisymmetric £ = 3-8 components, 
while the toroidal field is practically absent. 

Despite a consistent analysis using the same inversion 
methodology and MuSiCoS spectropolarimetric data of similar 
quality and phase coverage, the level of complexity turns out to 
be dramatically different for the surface magnetic field topolo- 
gies of a 2 CVn and 53 Cam. For a 2 CVn the I = 1 mode is 
more important relative to other components and the contribu- 
tion of the toroidal field is small. For 53 Cam the power is not 
concentrated in the dipolar component but spread out over the 
entire £ = 1-10 range, with the broad maximum at I — 4-7, and 
the toroidal components are noticeably stronger than in a 2 CVn. 
Therefore, we conclude that the surface magnetic field complex- 
ity of the two Ap stars studied using high-resolution four Stokes 
parameter observations is intrinsically different. 

6. Discussion 

In this paper we have described what is only the second self- 
consistent analysis of high-resolution Stokes IQUV line profiles 
of a magnetic Ap star. For the prototypical spectrum variable 
a 2 CVn, we have employed phase-resolved polarization spectra 
to derive a detailed map of the surface magnetic field intensity 
and orientation using the Magnetic Doppler Imaging approach. 
In addition, we have derived similarly detailed maps of the sur- 
face distributions of the abundances of Fe and Cr. 

We find that the overall structure of the magnetic field in 
a 2 CVn is dipole-like, with approximately half of the stellar sur- 
face covered with the outward-directed radial field while the 
other exhibits inward-directed field. This is in agreement with 
the field geometry suggested in previous magnetic analyses of 
a 2 CVn obtained from fitting of multipolar models. However, 
the field strength distribution we derive reveals that the field is, 
in fact, far more complex than the simple geometry suggested 
by earlier models. In particular, there is a definite asymmetry 
in the field strength and structure between the negative mag- 
netic pole and the positive pole. The field is clearly stronger at 
the positive pole and its structure is dominated by high-contrast 
magnetic spots where the field strength reaches 4.5 kG locally. 
Interestingly, this is approximately 1 kG higher than the local 
field at the opposite side of the star. These spot s are analogous 
to those detected in the A4p star 53 Cam by iKochukhov et al.l 
(2004b). We investigated the possibility that the defining char- 
acteristic of these maps - the high-contrast structure of the field 
strength at smaller spatial scales - may result from insufficiently 
smoothing as a consequence of underestimation of the regular- 
ization intensity. We thoroughly examined this possibility by 



conducting inversions with different intensities of the regular- 
ization parameter for the magnetic field. We confirmed that our 
choice of regularization is correct because any substantial in- 
crease of the smoothing of the magnetic maps noticeably wors- 
ened the agreement between the observed and computed lin- 
ear polarization observations. We thereby confirmed that the 
smaller-scale structures in our map are real, and required in or- 
der to reproduce the observations. We also computed magnetic 
field maps using just the Stokes IV spectra. These maps were no- 
ticeably less structured than those derived from the full Stokes 
IQUV data set, confirming that the linear polarization observa- 
tions are essential to detecting this complex, smaller-scale com- 
ponent of the magnetic field. 

A spherical harmonic decomposition was applied to quanti- 
tatively study the derived magnetic field topology. This analysis 
suggests that the field geometry of a 2 CVn is dominated by I — 1 
modes. At the same time, the contribution of the higher-/' modes 
and the amplitude of the toroidal field components is also non- 
negligible, with contributions of modes with t up to 5-6. 

We also re-visited the analysis of 53 Cam performed in 2004, 
obtaining a new field map of this star using the current version of 
the mapping software. Despite this consistent analysis using the 
same inversion methodology and spectropolarimetric data of the 
same origin, and of similar quality and phase coverage, the level 
of magnetic field complexity is found to be dramatically differ- 
ent for a 2 CVn and 53 Cam. For a 2 CVn the £ = 1 mode is more 
important relative to other components and the contribution of 
the toroidal field is small. For 53 Cam the power is not concen- 
trated in the dipolar component but is distributed over the entire 
£ = 1-10 range, with a broad maximum at £ — 4-7, and toroidal 
components that noticeably stronger than those of a 2 CVn. Also, 
the clear hemispheric asymmetry of the qualitative field proper- 
ties of a 2 CVn - with small-scale structures confined essentially 
to the stronger, positive pole, the origin of which is not known - 
is not obviously reflected in the field of 53 Cam. Therefore, we 
conclude that the surface magnetic field complexity of the two 
Ap stars studied using high-resolution four Stokes parameter ob- 
servations is intrinsically different. The reason for this difference 
is currently a mystery - similar maps of a larger sample of stars 
will be needed before a clearer understanding will be possible. 

The surface abundance maps of Cr and Fe show high- 
contrast distributions of both elements, varying in abundance by 
3-4 dex and with pronounced abundance minima at a location 
corresponding roughly to the negative magnetic pole. We also 
observed that the map obtained from the two strong Fe n lines 
exhibits noticeably smaller abundances in the region around the 
negative magnetic pole than those obtaine d for this star using 
weaker Fe lines by Kochukh ov et al.l d2002). We discussed how 
this could be be ascribed to the effects of vertical chemical strati- 
fication, which would lead to substantial weakening of the cores 
of intrinsically strong Fen lines relative to intrinsically weak 
lines, yielding a smaller abundance if such lines are analyzed 
neglecting chemical stratification. 

The investigations of th e magnetic fields of both a 2 CVn (in 
this paper) and 53 Cam ("bv IKochukhov et al.ll2004bl) lead to the 
view of a hidden complexity of the magnetic fields of Ap stars - 
a complexity which is only revealed with the help of linear polar- 
ization measurements. Despite the successes achieved, the data 
on which these studies have been based is fundamentally limited. 
The resolving power of the MuSiCoS spectropolarimeter was 
rather low - just 35000 - and the throughput of the instrument 
below 1%. These characteristics led to polarization spectra in 
which Zeeman signatures were often undetectable (particularly 
in linear polarization), or only detectable at relatively low signif- 
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icance. The current generation of high-resolution Stokes IQUV 
spectropolarimeters - including ESPaDOnS on the Canada- 
France-Hawaii Telescope, Narval on the Telescope Bernard Lyot 
(both with R * 6.5 x 10 4 and 20% throughput) and HARPSpol 
on the ESO La Silla 3.6m telescope (with/? * 1.2x1 5 and 10% 
throughput) - are capable of acquiring dat a of far greater qual- 
ity. In particular, observations acquired bv lSilvester et al.l {2008) 
with ESPaDOnS and NARVAL show the importance of high re- 
solving power for confidently measuring linear polarization. Not 
only are these instruments capable of acquiring much better data, 
they also allow us to acquire such data for a much larger sam- 
ple of stars. Application of these new tools to investigating the 
detailed magnetic topologies of Ap stars - with a range of field 
strengths, rotation rates, ages and masses - is critical to under- 
standing the physics of their magnetic fields. This includes not 
only the puzzling structures revealed in the present study, but 
also the much broader problems of field origin and evolution. 

In conclusion we would like to add that our studies of mag- 
netic field geometries of Ap stars unequivocally demonstrate the 
necessity of using spectra in all four Stokes parameters for re- 
liable reconstruction of the field topologies at various spatial 
scales. We found that a relatively simple and smooth field struc- 
ture inferred from circular polarization is invariably superseded 
by a substantially more complex magnetic topology when lin- 
ear polarization is incorporated in the Doppler imaging analysis. 
There are no reasons to believe this trend to be specific to Ap 
stars. Thus, we suggest that full Stokes vector observations of 
magnetic stars in other parts of the H-R diagram will lead to a 
similar dramatic change of the inferred field structure, especially 
for the late-type active stars which lack a dominant low-order 
field component. Understanding of stellar magnetism remains 
fundamentally incomplete without four Stokes parameter spec- 
tropolarimetry. 

Acknowledgements. Calculations presented in this paper were carried out at the 
supercomputer facility provided to the Uppsala Astronomical Observatory by 
the Knut and Alice Wallenberg Foundation and at the UPPMAX supercomputer 
centre at Uppsala University. O.K. is a Royal Swedish Academy of Sciences 
Research Fellow supported by grants from the Knut and Alice Wallenberg 
Foundation and from the Swedish Research Council. G.A.W.'s research is sup- 
ported by a Discovery Grant from the Natural Science and Engineering Research 
Council of Canada, as well as a grant from the Academic Research Program of 
the Department of National Defence (Canada). Resources provided by the elec- 
tronic databases (VALD, Simbad, NASA ADS) are gratefully acknowledged. 



Kochukhov, O., Bagnulo, S., Wade, G. A., et al. 2004a, A&A, 414, 613 
Kochukhov, O., Khan, S., & Shulyak, D. 2005, A&A, 433, 671 
Kochukhov, O. & Piskunov, N. 2002, A&A, 388, 868 

Kochukhov, O., Piskunov, N., Ilyin, I., Ilyina, S., & Tuominen, I. 2002, A&A, 
389, 420 

Kochukhov, O., Ryabchikova, T., & Piskunov, N. 2004b, A&A, 415, L13 
Kochukhov, O., Shulyak, D., & Ryabchikova, T. 2009, A&A, 499, 851 
Kochukhov, O., Tsymbal, V., Ryabchikova, T., Makaganyk, V., & Bagnulo, S. 

2006, A&A, 460, 831 
Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C, & Weiss, W. W. 

1999, A&AS, 138, 119 
Kurucz, R. L. 1993, SYNTHE spectrum synthesis programs and line 

data (Kurucz CD-ROM, Cambridge, MA: Smithsonian Astrophysical 

Observatory) 

Landolfi, M., Landi degl'Innocenti, E., Landi degl'Innocenti, M., & Leroy, J. L. 

1993, A&A, 272, 285 
Landstreet, J. D. 1982, ApJ, 258, 639 
Leroy, J. L. 1995, A&AS, 114, 79 
Lipski, L. & Stepieri, K. 2008, MNRAS, 385, 481 

Luftinger, T., Kochukhov, O., Ryabchikova, T., et al. 2010, A&A, 509, A71 
Maury, A. C. & Pickering, E. C. 1897, Annals of Harvard College Observatory, 
28, 1 

Mestel, L. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 

305, Astronomical Society of the Pacific Conference Series, ed. L. A. Balona, 

H. F. Henrichs, & R. Medupe, 3-15 
Piskunov, N. 2005, in EAS Publications Series, Vol. 17, EAS Publications 

Series, ed. G. Alecian, O. Richard, & S. Vauclair, 245-250 
Piskunov, N. & Kochukhov, O. 2002, A&A, 381, 736 
Pyper, D. M. 1969, ApJS, 18, 347 

Schaerer, D., Charbonnel, C, Meynet, G., Maeder, A., & Schaller, G. 1993, 
A&AS, 102, 339 

Schaller, G, Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 
Semel, M., Donati, J.-R, & Rees, D. E. 1993, A&A, 278, 231 
Shulyak, D., Tsymbal, V, Ryabchikova, T., Stiitz, C, & Weiss, W. W. 2004, 
A&A, 428, 993 

Silvester, J., Wade, G. A., Kochukhov, O, Landstreet, J. D., & Bagnulo, S. 2008, 

Contributions of the Astronomical Observatory Skalnate Pleso, 38, 341 
van Leeuwen, F. 2007, A&A, 474, 653 

Wade, G. A., Donati, J.-R, Landstreet, J. D., & Shorlin, S. L. S. 2000a, MNRAS, 
313, 851 

Wade, G. A., Donati, J.-R, Landstreet. J. D., & Shorlin, S. L. S. 2000b, MNRAS, 
313, 823 



References 

Babcock, H. W. & Burd, S. 1952, ApJ, 116, 8 

Bagnulo, S., Landi degl'Innocenti, M., & Landi degl'Innocenti, E. 1996, A&A, 
308,115 

Bagnulo, S., Wade, G. A., Donati, J.-R, et al. 2001, A&A, 369, 889 
Baudrand, J. & Bohm, T. 1992, A&A, 259, 711 

Bessell, M. 2000, Magnitude Scales and Photometric Systems, Encyclopedia of 

Astronomy and Astrophysics 
Borra, E. F. & Landstreet, J. D. 1977, ApJ, 212, 141 
Cohen, J. G. 1970. ApJ, 159, 473 

Cowley, A., Cowley, C, Jaschek, M., & Jaschek, C. 1969, AJ, 74, 375 

Donati, J.-R, Catala, C, Wade, G. A., et al. 1999, A&AS, 134, 149 

Donati, J.-R, Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, 

MNRAS, 291, 658 
Famsworth, G. 1932, ApJ, 75, 364 

Folsom, C. P., Wade, G. A., Kochukhov, O, et al. 2008, MNRAS, 391, 901 
Gerth, E., Glagolevskij, Y. V, Hildebrandt, G, Lehmann, H., & Scholz, G. 1999, 
A&A, 351, 133 

Glagolevskii, Y. V, Piskunov, N. E., & Khokhlova, V. L. 1985, Soviet 

Astronomy Letters, 11, 154 
Goncharskii, A. V, Ryabchikova, T. A., Stepanov, V. V, Khokhlova, V. L., & 

Yagola, A. G. 1983, Soviet Astronomy, 27, 49 
Khalack, V. R. & Wade, G. A. 2006, A&A, 450, 1 157 
Khokhlova, V. L. & Pavlova, V. M. 1984, Soviet Astronomy Letters, 10, 158 



12 



O. Kochukhov and G. A. Wade: Magnetic Doppler imaging of a 1 Canum Venaticorum in all four Stokes parameters 



Stokes I 



Stokes Q 



Stokes U 



Stokes V 



0.038 





^aTVvA, %A-^A^\/v 0706 




v^^^^v^ ^^^J^^\/^ °- 867 





0.949 
0.973 



I ■ I l.l.l L 



_L 



4824 4924 



J l.l.l L 
5018 4824 4924 



I I ■ I ■ I l.l.l L 
5018 4824 4924 



i.i.i i.i.i i.i.i i.i.i 
5018 4824 4924 5018 



Wavelength (A) 



Fig. 5. The comparison between observed (symbols) and synthetic (solid curves) four Stokes parameter spectra of a 2 C Vn. The 
thick curves show the best fit to the Fe and Cr lines obtained with the four Stokes parameter DI modeling of the magnetic field 
geometry (Fig. [6} and horizontal abundance distributions of both elements (Pig. [9]>. The thin lines show theoretical Stokes profiles 
illustrating the outcome of the MDI inversion with Tikhonov regularization enhanced by a factor of 10 relative to its optimal value 
(Fig. [8). Spectra for consecutive rotational phases are shifted in the vertical direction for display purposes. Rotational phases are 
indicated in the column to the right of the Stokes V panel. The bars at the lower left of each panel show the horizontal and vertical 
scale (0.5 A and 1% of the Stokes / continuum intensity respectively). The number in brackets at the top of each panel quotes the 
mean deviation between observations and spectrum synthesis for that Stokes parameter. 
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Fig. 6. Surface magnetic field distributions of a^CVn derived from the Stokes IQUV profiles of the Fen and Crn lines. The star 
is shown at five equidistant rotational phases as indicated at the top of the figure. The aspect corresponds to the inclination angle 
i = 120° and vertically oriented rotational axis. The upper row of spherical plots visualizes the distribution of field strength, with 
contours of equal magnetic field strength are plotted every 0.5 kG. The thick line shows the stellar rotational equator. Rotational 
axis is indicated with the vertical bar. The lower panel shows the orientation of the magnetic vectors. In these vector maps the black 
arrows show field vectors pointing outside the stellar surface and the light arrows correspond to the vectors pointing inwards. The 
arrow length is proportional to the field strength. 
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Fig. 7. Rectangular projection of the magnetic field geometry of a 2 CVn derived from the Stokes IQUV profiles of the Fe n and Cr n 
lines. The longitude coordinate divided by 360° corresponds to the phase of the subsolar meridian in Fig. [5] The dotted line shows 
the highest visible latitude for the adopted inclination angle i - 120°. The lines of equal field strength are shown for the [-3, +3] kG 
interval, with a 1 kG step. 
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Fig. 8. Same as Fig.[6]but for the Stokes IQUV imaging with ten times larger Tikhonov regularization for the magnetic field. 
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Fig. 9. Surface distribution of the Cr and Fe abundances derived simultaneously with the magnetic field geometry in Fig. [6] The star 
is shown at five equidistant rotational phases as indicated at the top of the figure. The aspect corresponds to the inclination angle 
i = 120° and vertically oriented rotational axis. The element concentrations are given in the logarithmic units relative to the total 
atomic number density: e(El) = log (A^i/MoO-The contour lines in the spherical abundance maps are plotted with a step of 0.5 dex. 
The thick line shows the stellar rotational equator. Rotational axis is indicated with the vertical bar. 
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Fig. 10. Same as Fig.[5]except that theoretical four Stokes parameter spectra shown with thin lines correspond to the best-fit magnetic 
geometry reconstructed in the Stokes IV imaging with multipolar regularization. The corresponding magnetic field geometry of 
a 2 CVn is presented in Fig.fTTI 
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Fig. 11. Same as Fig.[6]but for the Stokes IV imaging with multipolar regularization. 
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Fig. 12. Coefficients for the spherical harmonic expansion of the MDI maps of a 2 CVn and 53 Cam. The columns correspond to the 
surface field of a 2 CVn inferred from the full Stokes vect or inversion (left), magn etic field geometry of a 2 CVn obtained from the 
Stokes IV spectra (middle) and magnetic field of 53 Cam dKochukhov et al. 2004a). The top row shows coefficients for the poloidal 
field, while the bottom row corresponds to the toroidal field component (practically absent in the MDI maps of a 2 CVn reconstructed 
from the circular polarization spectra). 



